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Abstract 


EXPERIMENTAL INVESTIGATION OF A 
2 1/8"-DIAMETER CONSTANT-AREA AEROTHERMOPRESSOR 
WITH SUPERSONIC INLET 


by 
Robert Torrey MacKay 


Submitted to the Department of Mechanical Engineering on 
May 27, 1955 in partial fulfillment of the requirements for 
the degree of Master of Science in Mechanical Engineering. 


It is known from theoretical studies that good per- 
formance of the Aerothermopressor, currently under study 
at the Massachusetts Institute of Technology under the di- 
rection of Professor A. H. Shapiro and under the sponsor- 
ship of the Office of Naval Research, is dependent upon the 
ability to produce extremely small droplets upon injecting 
cooling water into a hot, high-velocity gas stream. 


The Aerothermopressor, a circular duct fitted with 
a water injection nozzle, produces a rise in the stagnation 
pressure of such a gas stream by the mechanism of evapora- 
tive cooling, and has as its main objective the improve- 
ment of gas turbine plant performance. To accomplish this 
objective it-must necessarily rely upon a high evaporation 
rate. 


By taking in gas-turbine exhaust gases, raising their 
stagnation pressure, and exhausting to atmospheric pressure, 
the Aerothermopressor will in effect produce a vacuum in 
the turbine exhaust region, with consequent increase in tur- 
bine power and efficiency. 


The need for small droplets with high surface-volume 
ratio for higher heat transfer and evaporation rates has 
given rise to an additional research program at the Massa- 
chusetts Institute of Technology for the purpose of devi- 
Sing a suitable method for accurately measuring the drop 
sizes of water droplets moving in a high-velocity gas strean, 
and ultimately to determine by drop-size measurements what 
type of water injection nozzle and what dynamic and thermo- 
dynamic conditions will produce the smallest droplets. 


At the present time the belief is being adhered to- 
that droplet size under conditions existing in the Aero- 
thermopressor is sensitive to the difference between inlet— 
gas velocity and water injection velocity, the higher rela- 
tive velocity producing the smaller droplets. 


LE 770 


Aerothermopressor evaporation takes place as the wa- 
ter droplets travel downstream with the gas, with accompany- 
ing decrease in the temperature gradient between the two. 


With an upper limit on turbine exhaust (Aerothermopressor 
inlet) stagnation temperature being imposed by turbine power re- 
quirements, the steady flow energy equation suggests that 
with both high gas velocity and high gas stream temperature 
being favorable to high evaporation rates, an optimum inlet 
velocity must exist for a given inlet stagnation temperature 
and water temperature. 


It was the purpose of this thesis to determine experi- 
mentally whether this optimum velocity lies above or below 
the speed of sound for a typical turbine exhaust stagnation 
temperature of 1500°R. 


Tests of a small-scale Aerothermopressor were conducted 
at inlet Mach numbers of approximately 1.5 and 1.35. The 
results were compared with subsonic data previously recorded 
for the same rig. From this comparison the conclusion was 
drawn that the optimum inlet velocity for the Aerothermopressor 
operating at 1500°R lies in the supersonic range. 


The main design feature incorporated in this thesis was 
to devise a means of converting a bellmouth nozzle to a super- 
sonic nozzle. This was accomplished by the design of a water 
injection nozzle se shaped as to also function as an area plug 
to give an annular-shaped converging diverging nozzle. 


Thesis Supervisor: Ascher H. Shapiro 
Title: Professor of Mechanical Engineering 
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I. ginGreduction 

It is known (Reference 1) that one of the effects of 
cooling a gas stream is to raise its stagnation oressure. The 
"NerothermoPpressor (Reference 2) was conceived by Professor 
A. H. Shapiro at the Massachusetts Institute of Technology 
to exploit this effect by injecting cooling water into a hot 
Pe weuoci ty gas stream. With the advent of the gas turbine 
as a Prime mover for ship propulsion and power generation, 
interest has focused upon the Aerothermopressor (abbreviated 
"A.T.P!) as a means of improving gas turbine performance. 

By receiving gas turbine exhaust gases, injecting 
cooling water into them, and exhausting to atmosphere the 
Aerothermopressor, with its attendant rise in stagnation 
pressure from inlet to exhaust, will provide the turbine 
with a back-pressure below atmospheric. Thus the turbine 
exhausts into a vacuum, with a consequent increase in its 
power and efficiency. The Aerothermopressor can be said to 
perform for the gas turbine cycle the same function as the 
cOndenser does for the steam turbine cycle, the only operating 
cost being the small amount of power required to pum sea 
water or river water into the Aerothermopressor at low velocity. 

It appears from present knowledge of the phenomena 
occurring in the Aerothermopressor process that the best per- 
formance will eventually be achieved with either a high sub— 
sonic or low supersonic gas inlet velocity. The effect of 


inlet Mach number upon A.T.P. performance is discussed in 


—_— st 





section Lil. 

The term "Supersonic inlet" is utilized rather than 
the term "supersonic Aerothermopressor" since it has been 
found (Ref. 2) that in the. subsonic-inlet A.T.P. process a 
continuous transition from subsonic to, supersonic‘ velocity 
is caused by evaporation effects. 

The experimental and theoretical research program 
now in progress at the Massachusetts Institute of Technology 
under the sponsorship of the Office of Naval Research has 
thus far concentrated upon subsonic gas inlet velocity, with 
its advantage of greater experimental flexibility and economy 
(variation of gas flow and inlet Mach number in the supersonic 
range ordinarily requires either a sevarate nozzle for each 
Mach number, or flexible walls in the nozzle region). 

This thesis was undertaken for the purpose of ob- 
taining some experimental data at "low" supersonic inlet 
velocities for comparison with data already recorded at high 
subsonic inlet velocities, to aid in answering the question 
as to whether optimum verformance will demand an inlet gas 
velocity above or below the speed of sound. The only super- 
sonic inlet A.T.P. data reported at the time of this writing 
was obtained at approximately Mach 2.0 on a 1.525" diameter 
constant-area A.T.P. without diffuser (Ref. 3). 


At the time of this investigation a "medium" (11" dia- 





meter test section) constant-area A.T.P. in the Gas Turbine 
Laboratory at the Massachusetts Institute of Technology was 
being used for subsonic-inlet experimental investigations of 
the effect of area variation upon performance, by means of 
conical area-plugs inserted inside the test section and moved 
longitudinally along the axis. This means of area-variation 
was chosen (with the disadvantage of added stagnation-pres- 
sure loss due to drag at the surface of the area plugs) over 
the more desirable possibility of varying the test section 
Giameter along the length, due to greater flexibility and 
economy. 

Prior to the development and installation of this 
"medium-scale" rig, Wadleigh in his doctoral dissertation 
(Ref. 4) investigated subsonic-inlet A.T.P. performance in 


' diameter constant- 


a "small-scale" rig consisting of a 2 1/ 
area test section supplied with hot gas from standard 6" pipe 
through an elliptical-contour bellmouth nozzle, and exhausting 
through a constant-angle diffuser into 6" standard pipe. 

Since Wadleigh's apparatus was available for use, and since 
the very nature of the Aerothermopressor process (high temp- 
erature plus corrosive gases) calls for stainless steel con- 
struction throughout, it was deemed advantageous from the 
point of view of economy as well as time to utilize this 
apparatus in the supersonic-inlet investigations. 


At small diameters the predominance of wall friction 


effects (tending to reduce the stagnation pressure) preclude 








the possibility of obtaining an actual rise in stagnation 
pressure by evaporative cooling, but as Wadleigh has pointed 
out, the measure of amount of reduction in stagnation vressure 
loss obtainable by introducing cooling in a small-scale rig 
may be used with some success as a criterion for estimating 
A.T.P. verformance. 

A complete summary of the exnerimental and theoretical 
work accomnlished on the Aerothermovressor project to date 
(as well as a complete bibliography on the subject) is con— 


tained in Reference 2. 
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II. Main Test Equipment and Measurements 

A schematic drawing of the main test equipment used 

in this investigation is given in Fig. 1. For reference 
purposes, all pressure taps are numbered in order, commencing 
at the upstream end of the apparatus. Two possible combi- 
nations were afforded by Wadleigh's apparatus: a constant- 
areatest section 72" long in conjunction with a 6-degree 
diffuser, or a constant-area test section 36" long in con- 
junction with a 3-degree (total included angle) diffuser. 
Due to the problem of friction choking in a constant-area 
test section at low supersonic inlet Mach numbers, the 
shorter test section (2.125" I.D., 36" long) with the 
3-degree diffuser (approximately 78" long) was chosen for 
the supersonic-inlet testing. 

The test section was made from heavy-wall #321 stainless 
steel tubing and fitted with eleven static pressure taps 
.030" in diameter, spaced at intervals along the length 
(taps #2 to #12 in Fig. 1). The diffuser was rolled from 
#321 stainless steel sheet 1/8" thick and fitted with ten 
static pressure taps .030" in diameter spaced at intervals 
along the length (taps #13 to #22 in Fig. 1). The ellip- 
tical bellmouth was machined in a #321 stainless steel block 
1 1/2" thick, welded to the inlet end of the test section. 
Table I tabulates the axial locations of all pressure taps 


and other pertinent stations along the duct length. 
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; Hot combustion gases were supplied to the bellmouth 
through 6" stainless steel pipe by a natural gas /"city gas") 
furnace constructed by a the Etter Engineering Co., and 
utilizing an Eclipse NHE burner No. 5 as the primary heat- 
ing burner and an Eclipse Walltite LEA 9 burner No. 3 as 
the pilot burner. Air was supplied to the furnace by a 
centrifugal blower. The gas-air mixture to the furnace was 
fixed by a pressure-regulator control system: hence good 
temperature control was achieved by merely controlling the 
air supply to the furnace with a large butterfly valve 
and a small bypass valve located in the air supply line 
between the blower and the pressure regulator. 

The initial stagnation temerature Ty), of the gas 
was measured inside the standard 6" pipe at a location 9" 
from the entrance of the bellmouth, by means of a five 
shielded chromel-—alumel thermocouple manufactured by the 
Aerotech Specialties Co. in Glastonbury, Conn., and read 
by a Leeds and Northrup K-2 potentiometer. 

The initial stagnation pressure po, was measured by 
means of a stagnation pressure probe (tap #1 in Fig. 1) lo- 
cated in the 6" hot gas supply line 15" upstream from the 
bellmouth entrance, 

Before exhausting the hot gases into the campus 
atmosphere, further cooling not necessary to the A.T.P. 


process was effected by discharging the diffuser through 
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a two-foot length of standard 6" pipe into a quench tank 
7' high and 3' in diameter, fitted inside with two water 
sprays. Final stagnation pressure Pop was measured by 
means of a stagnation vrobe (tap #23 in Fig. 1) located 
in the 6" discharge pipe. To obtain and maintain gas 
flow through the test section, suction was provided at 
the top of the quench tank by means of a 6" steam ejector 
which is vermanently installed in the Gas Turbine Labo- 
ratory. Back vressure on the diffuser exit was controlled 
by a 6" gate valve located between the 6" steam ejector 
and the quench tank. Excess quench water was pumped out 
of the bottom of the quench tank by adapting a small in- 
jector (which was readily available) to the job of an 
ejector. 

All pressures were measured on a mercury mano- 
meter board as differences (cm. Hg.) from an atomspheric 
mercury column, and converted to absolute pressures by 
the local mercury barometer reading. 

The problems of (1) furnishing the A.T.?. cooling-— 
water supply at the entrance to the test section and (2) 
adapting Wadleigh's bellmouth nozzle to supersonic flow 
were solved simultaneously by designing one niece of equin- 
ment which for lack of a better name was called a "water 
injection nozzle assembly". Details of this assembly are 
ghown in Figs. 2 and 3. 


The main body of the water injection nozzle assembly 
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functions as an area plug which, when mounted inside the 
bellmouth and test section, yields an annular-shaped con- 
verging-diverging supersonic nozzle. The area variation 

of this plug was designed on a One-dimensional isentropic 
basis so as to fix the throat of the annular-shaped suner- 
sonic nozzle at the exit plane of the bellmouth nozzle. since 
a static pressure tap (tap #2) was located at this point. 

In keeping with the current line of thinking in the 
Aerothermovressor Project as to means of water injection, 
the water injection nogzle supplies its water through six 
"atomizer tubes" of standard 12-gauge stainless steel tubing 
manufactured by the Hub Needle Company, Boston, Mass.. dischar- 
ging parallel to the direction of gas flow. These six atom — 
izer tubes were spaced symmetrically about the test section 
axis at the root-mean-square radius of the test section so 
as to serve equal gas flow areas. Fig. 4 shows the water 
injection nozzle assembly mounted in place inside the bell 
mouth and test section to form the annular-shaped sunersonic 
nozzle. 

With this scheme, variation in inlet Mach number is 
obtainable by merely changing conical tins in the nozzle 
assembly. Two tips were actually manufactured and tested: 
conical tip No. 2 (Fig. 5) with diameter .840" at bell- 
mouth exit plane (designed for Mach 1.5), and conical tip 


No. 3 (Fig. 6), with plug diameter .540" at bellmouth exit 
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plane (designed for Mach 1.35). Tip No. 1 was never manu- 
factured, having been abandoned in favor of the smoother 
contours of tip No. e. 

The water injection nozzle assembly was manufactured 
of #304 stainless steel by the instrument room of Shop #31 
at the U.S. Naval Shipyard, Boston, Mass. 

Since the conical tip (and therefore the nozzle sec- 
tion) extends 4" into the 36" test section, the actual ef- 
fective length of the constant-area test section for the 
supersonic runs is 32". 

Water was supplied to the water injection nozzle 
assembly at city water pressure, and water flow was measure 
by a Schutte and Koerting rotameter (rotor No. 4B, tube 
No. 4R). 

Gas flow through the choked converging-diverging 
nozzle was estimated by "Fliegner's formula" (with k = 1.35), 
modified by an estimated nozzle discharge coefficient 
C, - -98. For this purpose the furnace combustion gases 


W 
were assumed to have the same molecular weight and specific 
heatsas air. The ratio of speeifice heats k= 1.35 was 

chosen as a mean value for the range of stream temperatures 
encountered in the two nozzles. For ready reference, isen- 


trepic flew Bunetions for k = 1.35 are tabulated in Table II. 

















‘* ; 
—<=— @! €faI Gees ae @ 
- 
- A ieee 8 Pome a a 
(9 om SO eee ol Qo Gee 
+ — —<—s TT 
-—ait = Oe ae! Lorem al 
> a = oa Ge G at 
ee eet ee ee 
— 
—=- - —-— | == ow 
——_ eT ee —— — | a al ie 
_» se 4 — bt =—t oe eee 












a => 














<a _ .——i ale oe 

“ re —? @ Jae te 6 
- - ~_i-_- & «& Le 

“a ioeeke = = @& 

















* 
= - e = * aay 


———— ei ee ew SS SS 


LIL. Effects of Inlet Mach Number on A.T.P. Performance 

Based on a one-dimensional analysis utilizing the 
continuity, momentum, and energy equations for a perfect 
gas and neglecting the less influential effects of: (1) 
change in mass rate of gas flow due to "new water vavor 
entering the gaseous phase by evaporation from the water 
droplets, (2) change in molecular weight of the gas due 
to the appearance of this same new water vapor, and (3) 
change in the ratio of specific heats k due to this new 
vapor and due to temperature change, the governing dif- 
ferential equation for the local stagnation pressure po 
of the gas at any cross-sectional plane through an A.T.P. 
duct without shock may be written (Ref. 2) : 


dp, =~ Apa , 4¢ 22 +n) 


lo 


In this equation M and cr represent the local Mach number 
and stagnation temperature, z the axial distance along the 
duct from any Fixed reference point (z increasing positive 
in the direction of the flow), D the local diameter, f the 
Fanning friction factor, and n a term representing drovlet 
drag exerted on the gaseous phase. A more refined equation 
which includes the effects neglected here (or, for that 
matter, the governing differential equation for any of the 
other dependent flow variables) may be written down at 


once by referring to the Table of Influence Coefficients 
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in Ref. 2; but the above simple form is adequate for the 
purpose intended here. The above equation will be referred 
to henceforth as the fundamental governing equation of the 
ao. P. 

The factor enclosed in parentheses on the right-hand 
side of the fundamental equation represents the net influence 
of the major effects in the A.T.P. process: (1) evanorative 
cooling (decrease in gas a ), tending to increase the stag- 
nation pressure, and (2) wall friction(assisted by the 
initial water droplet drag), tending to decrease the stag- 
nation pressure. In actual full size A.T.P. operation in 
which evaporative cooling effects vredominate, the in- 
tegrated form of the factor in parentheses will be nositive, 
yielding a net rise in stagnation pressure from inlet to 
exhaust. Further, it may be seen from the factor me that 
for such full-size overation the higher the Mach number 
level through the duct, the higher will be the stagnation 
pressure rise. For supersonic oneration with shock, the 
loss in stagnation pressure across the shock must also be 
reckoned. In constant-area operation this shock influence 
must necessarily limit supersonic operation to the lower 
end of the supersonic range: with variable area in the 
duct (studied theoretically by Gavril, Ref. 7, and currently 
under experimental study with subsonic inlet by Assistant 


Professor A. A. Fowle) the possibility exists of diffusing 
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supersonic flow to a lower Mach number before shock, thus 
reducing shock losses. 

The best A.T.P. performance does not, however, lie 
with the highest possible inlet Mach number as it might ap- 
pear from this reliminary discussion, since the factor in 
parentheses in the above equation is also a function of Mach 
number. The term dT, /T, is the change in gas stagnation 
temperature brought about by evaporation of the water drop- 
lets, which in turn is greatly influenced by inlet Mach 
number, as discussed in the following paragraphs. 

Water is believed to leave the "atomizer tubes" of 
the water injection nozzle in the form of "ligaments" or 
"sheets" which are broken up by the impact or drag of the 
high-velocity gas whto water droplets (Refs. 5 and 7). It 
is further believed that the size of water droplets re- 
sulting from this process depends in large measure upon the 
relative velocity between the gas the the water--the higher 
the relative velocity the smaller the drops (Ref. 7). So 
important is drop size to A.T.P. performance that an addit- 
ional research program is underway at the Massachusetts In- 
stitute of Technology on droplet technology, intended first 
to investigate several schemes of measuring drop size in 
a moving stream, and eventually to determine what conditions 
will yield the smallest drops. Information obtained from this 
program should be of value to combustion studies as well 
as to A.T.P. design. 


For a fixed water-to-air ratio and a fixed initial 





temperature difference between water drovlets and gas, the 
smaller the droplet diameter the greater is the amount of 
water surface area exposed for heat transfer (by simple 
geometrical considerations) and hence the greater the 
evavoration rate. In addition to this effect of increased 
heat-transfer area, a smaller diameter means a lower Rey- 
nolds number: hence a higher coefficient of heat transfer. 
Gavril (Ref. 7) has demonstrated that due to the combi- 
nation of these two effects the heat transfer rate between 
water droplets and gas varies inversely with roughly the 
square of the droplet diameter. Since evaporation due to 
heat transfer is the heart of the A.T.P. vrocess. the im-— 
portance of obtaining small water-dronlet diameters can 
not be overemnhasized. 

The evaporation must also devend unon the difference 
between gas stream temperature and water droplet temnerature. 
During the evaporation process the water droplets are heated 
un and the gas is cooled: hence the temnerature difference 
between the two is gradually reduced as the moist gas and 
water drovlets move downstream. To maintain a high evano 
ration rate, then, a large initial temperature difference 
between air and water (at the water injection point) is 
desirable. 

For a fixed initial stagnationtempverature (such as 
the design turbine exhaust temperature of a gas turbine de- 


signed to be fitted with an Aerothermopressor) the gas may 
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be accelerated adiabatically to a high subsonic or low 
supersonic velocity prior to entering the A.T.P.; but (by 
the energy equation for steady flow) as the velocity is 
increased the stream temperature of the air is decreased. 
Thus the ontimum inlet Mach number for the A.T.P. must 
achieve a balance between two contradictory requirements: 
high gas velocity to obtain small droplets from the water 
"preakup" process, and high gas stream temperature, both 


of which are favorable to high evaporation rates. 
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IV. Limitations Imposed by Test Section Length 

Constant-area Aerothermopressor operation with 
supersonic inlet implies the existence of a shock some- 
where in the duct. Since the purpose of the A.T.?. is to 
obtain a rise in stagnation pressure and since shock from 
high supersonic Mach number produces large losses in stag- 
nation pressure, good supersonic-inlet A.T.P. operation 
at constant area must necessarily be limited to low suver- 
sonic Mach numbers with their attendant "weak" shocks. 

For this reason it was considered desirable to in- 
vestigate suversonic Mach numbers of about 1.5 and below. 
How ever, the length of the test section (length L = 32", 
diameter D = 2.125", L/D = 15) imposed a minimum value on 
the supersonic Mach numbers which could be achieved. 

For a given constant~—area test section with Fanning 
friction f and a length-diameter ratio L/D, there is in 
"dry" gas flow a maximum subsonic and minimum supersonic 
Mach number which can exist at the test section entrance, 
corresponding to friction choke at the test section exit. 
(Ref. 8) These values are tabulated in the "Fanno Line" 
tables (eg. ref. 6). With the 15-diameter length of test 
section employed in these tests assuming, for instance, a 
friction factor f - .005, the computed value 4fL/D = 0.3 
yields (for k - 1.4) values of Msubsonic_., = 0.66 and 


1.98. 


M ‘ = 
supersonicisy 





ae 


Thus it was not to be exvected that dry sunersonic 
flow could be achieved in the existing test -section at Mach 
numbers 1.50 and 1.35, and dry supersonic flow was in fact 
found experimentally to be impossible. (Fig. 9). 

The theoretical one-dimensional analyses which have 
been carried out (Refs.4 and 7) have resulted in a "Table 
of Influence Coefficients" (see also Ref. 2) which show the 
effect of a change in any one of the arbitrarily chosen 
independent properties upon each of the remaining devendent 
properties associated with A.T.P. flow. Use of these co- 
efficients in interpreting Wadleigh's experimental data 
and Gavril's numerical computations (carried out on the 
"Whirlwind" high-sveed electronic digital computer at 
Massachusetts Institute of Technology) has resulted in 
Table III (reproduced from Ref. 2) which shows the effects 
of the four major controlling parameters of the Aerothermo- 
pressor (area change, evaporation, wall friction, and 
liquid acceleration) upon the flow properties. 

From Table III it may be seen that the effect of 
wall friction in the A.T.P. is to "drive" the Mach number 
toward unity (just as in dry Fanno-type flow), while the 
effect of evaporation is to drive the Mach number away 
from unity. Knowing these two effects to be opposite to 
each other, it was expected (and actually realized exveri- 
mentally) that "wet" A.T.P. supersonic flow might be ob- 


tained at supersonic Mach numbers lower than the limiting 
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value imposed by friction choke in the dry flow. Conical 

tip No. 2 (Mach 1.5) was tried first and found to provide 

wet supersonic flow in the test section. After this was 
found successful, a still lower Mach number of 1.35 (coni 

cal tip No. 3) was tried and found also to be successful. 

Mach number 1.35 is believed to be very close to the minimum 
supersonic inlet Mach number for which supersonic-inlet A.T.P. 
operation can be achieved with the present test section since 
at this Mach number supersonic flow was obtainable only at 
avery critical value of water air ratio. Any slight in 
crease or decrease in the water flow from this level resul- 
ted in unstable flow which quickly resulted in steady sub- 
sonic flow with friction choke at the test section exit. No 
attempt was therefore made to go to supersonic Mach numbers 
below 1.35. 

In his subsonic-inlet runs Wadleigh was able to esti- 
mate nozzle losses and Fanning friction factor by "hot dry" 
and "cold dry" data. This was not possible for the superson- 
ic-inlet case. since dry supersonic runs were prohibited by 
the friction-choke phenomenon just mentioned. However, at 
Wadleigh's suggestion, with the water supply to the Water 
Injection Nogzle cut off, auxiliary water was injected ra- 
dially into the stream through six symmetrically-spaced 
taps in the test-section wall located 6" downstream from 


the nozzle exit (midway between pressure taps #5 and #6). 
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The evaporative effects of this auxiliary water were 
sufficient in the case of the Mach 1.5 nozzle to eliminate 
the friction choke at the test section exit and produce 
supersonic flow in the entire test section. This gave 
"dry" supersonic pressure data as far downstream as tap #5 
(Run #34, Fig. 9), from which actual flow Mach No. and 
nozzle losses were estimated through the well-known for- 
mulas of "dry" gas dynamics (Section VI). In the case of 
the Mach 1.35 nozzle, this radial water injection was not 
sufficient to overcome the friction choke; hence dry super- 


sonic data for this nozzle is not available. 
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V. Experimental Results 

For the purpose of comaring Aerothermopressor per- 
formance at low iseskdtli’c entry with performance at sub- 
sonic entry, Aerothermopressor runs (runs with water in- 
jection at test-section inlet) were made at nominal inlet 
Mach numbers 1.5 and 1.35, utilizing conical tins No. 2 
and No. 3. An initial stagnation temmerature TO, = 1500°R 
was chosen for this comparison for two reasons: (1) it 
represents a reasonable exhaust temperature for an open- 
cycle gas turbine and (2) considerable subsonic-entry data 
was available at this temperature for comparison purposes. 

Curves of static pressure vs. axial distance along 
the duct for these runs are plotted in Figs. 7 and 8. Since 
the velocity at the stagnation pressure taps #1 and #23 
(inside the 6" pipe) is negligible, the stagnation pressures 
measured at these two stations also represent (within the 
accuracy of the mearsurements themselves) the static pvres- 
sures. Since the walls of the furnace were partially open 
to atmosphere,the initial stagnation pressure Po, at tap #1 
was essentially atmospheric in all the runs. For conven- 
ience, the pressuresare plotted in the dimensionless ratio 
form D/P; » and a schematic of the duct (showing locations 
of the various pressure taps by reference number) is drawn 


directly on the graphs. 


Fig. 9 presents the dry characteristics (no water 
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supplied to the stream) of the apparatus with the two dif- 
ferent nozzles, showing for the Mach 1.5 tip a shock in the 
supersonic nozzle followed by subsonic friction choke at 

the test-section exit, and for the Mach 1.35 nozzle a purely 
subsonic flow with friction choke at the test-section exit. 
As mentioned previously, Fig. 9 also shows the dry super- 
sonic nozzle pressure data obtained for the Mach 1.5 nozzle 
by downstream auxiliary radial injection. 

As mentioned in Section IV, the water-air ratio re- 
quired to maintain wet supersonic flow with the Mach 1.35 
nozzle was extremely critical; hence water-air ratio was 
not available as a variable. With the Mach 1.5 nozzle, the 
band of water-air ratios for which stable supersonic flow 
could be maintained was slightly wider. With this nozzle 
runs were made at approximately the maximum and minimum 
water-air ratios permissible for gupersonic flow, but even 
here the band was still so narrow that vractically no 
effect on performance could be imposed by varying the water 
suoply.. 

For the runs in Figs. 7 and 8 the water supply valve 
was adjusted to a point midway between the two cut-off 
points (too much water and too little water) at which the 
supersonic flow became unstable. The water-air ratios for 
this: condition Se Sowa vet were .235 for Mach 1.5 and .169 
for Mach 1.35, which may be considered to be "optimum" value 


for the particular apparatus used. 
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Before 1500°R was selected for the supersonic-subsonic 
comparisons, the initial stagnation temperature Or was 
varied between the minimum temperature required for stable 
operation and a maximum of 1800°R. Within this range of 
temperatures the effect of To; on A.T.P. verformance was 
very slight-~-the variation in pressure differences obtain-—- 
able by varying Tos being measured in fractions of a centi- 
meter of mercury. 

in intervreting thesolo@seof Figswe7,S, and 9 it 
should be borne in mind that the curves were constructed 
by simnly connecting with straight-line segments the 
pressure values which were reasured at intervals along 
the duct length. For this reason the peaks and valleys 
in these plots do not necessarily represent the maximum 
and minimum pressures in the actual flow--these may just as 


likely have occurred between vressure tans. 
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VI. Determination of Inlet Mach Number and Nozzle Losses 
Although conical tips No. 2 and No. 3 are referred 
to throughout this report as providing "Mach 1.5" and 
"Mach 1.35", respectively, these are merely nominal Mach 
number values. The geometric area ratios for the annular 


nozzles formed by these tips are representable in a one- 


dimensionable flow analysis by: 


AK 
A sest a Dest 
oe = = s 
A throat D yest a D tip. rit ori 6 Dee 


where Digg; is the diameter of the test section (2.125"), 
Dtipopit ts the "critical" tip diameter at the cross-section 
lying in the plane of the bellmouth exit (.840" for tip No. 
2: .540" for tip No. 3), and Diypeis the outside diameter of 
the atomizer tubing (.109"). Since the diameter of the 
Water Supply Plug (Fig. 2, part A) is .840", the annular 
nozzles formed by tips No. 2 and No. 3 are of the character 
illustrated in the sketches below, where the elliptical sr 
tour has been replaced by a circular contour and the tips 
are represented by large-angle right-circular cones to exag. 
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From these sketches it can be seen that while the 
geometric throat formed by tip No 2 is theoretically lo- 
cated exactly at the bellmouth exit, the throat formed by tip 
No 3 is actually slightly upstream of the bellmouth exit; 
hence the foregoing revresentation for nozzle area ratio 
is somewhat less securmate for tip No.3 than for tip No 2. 
This displacement is probably not great enough to intro— 
duce any serious error in the computed value of area ra- 
tio in this nozzle; but it must be borne in mind in inter- 
preting the curves of Fig. 8 that tap #2 does not record 
the vheget pressure for this nozzle. This accounts for 
the peculiar characteristic of the pressure curves between 
taps #2 and #3 in Fig. 8. 

Based upon this method of computing area ratio. the 
area data and the corresponding values of Mach number M 
from the isentropic flow functions of Table II are tabula- 
ted below for the two nogzles. The area of the test section 


is 3.547 squares inches. 


Nozzle Nozzle 
Conical Nominal Throat Area Area Corresponding 
Tip Mach No. (sq. in.) Ratio M (isentropic) 
No. 2 ee 2.936 1.208 1,35 
NO. 3 Ls) 3.263 1.087 1.34 


With the dry nozzle pressure data obtained for the 
Mach 1.5 tip by auxiliary radial water injection downstream 
of the nozzle (Section IV and Fig. 9), a very convenient 


method (Ref. 8) is available for calculating the nozzle exit 
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Mach number, utilizing an estimated nozzle discharge coeffici- 


ent Cy: | | 
(BAN. =e Clas, EE geom 


In this relationship the quantity on the left repre- 
sents a commonly tabulated (e.g., Ref. 6) isentronic flow 
function. For easy reference, this function is tabulated 
in Table II for k - 1.35. With the measured nozzle exit 
- .274 (tan #4, run #34) and the 


pressure ratio D/Do4 


above-recorded geometric area ratio: 


(Olle = TANI. 3g 


From Table II this gives a nozzle exit (test-section 
entrance) Mach number of 1.47 for conical tip No. 2. 

Taking this value of 1.47 to be the best available 
estimate of the actual flow Mach No. at nozzle exit, an 
estimate may be made of the nozzle losses. Denoting 


nozzle exit conditions by the subscript "e": 


fea= lo] gives Pe/Po, = .290 (From Table IT) 
De/Po, = .274 (measured, tap #4. run #34) 
De/Po 
Pog/Po, = ———— = .274/.290 = .945 
Pe Po, 


It should be noted that this loss parameter is based 


upon an assumed value of .98 for the nozzle discharge coef— 





ficient. 
The rise in stagnation pressure exhibited in the dry 
constant-area supersonic flow between tans #4 and #5 in run 


#34 is attributed to wall friction (Table IIT). 
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VII. Mass Flow of Air 

Following the terminology which has come into usage 
in the Aerothermopressor Project, the hot dry furnace gases 
entering the A.T.P. will be referred to henceforth in this 


tt 


report as "air", the term "gas" being reserved for the gase- 
ous phase of the wet mixture of combustion gases, water va- 
por, and water droplets which exists in the duct after water 
injection. | 

Mass flow of air was calculated from Fliegner's for- 


mula (Ref. 8) for a choked nozzle, modified by a nozzle dis- 


charge coefficient Cw: 









ke 
SS (fe) 


Taking k = 1.35, R= 53.35 ft-lbf , go= 32.17 lbm2ft 
Tbm-OR Tbf-sec 





Po; Athroat 


Vie 


and C, = .98, this has the form: 


Wa = . 515 Po: A throat 
mt, 
At the measured operating values of oor 14.5 psia 
and be 1500°R, this gives for the air flow passed by the 
two nozzles: 


Air Mass Flow We 
Conical Tip (1b. /sec. ) 


No. 2 tae lg) 565 
No. 3 (Mach 1.35) .629 





VIII. Interpretation of the Static Pressure Curves 
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The three major regimes (Ref. 2) of constant-area 
A.T.P. behavior are clearly visible in run #32 (Fig. 7) 
and run #44 (Fig. 8). 

In regime I droplet drag is predominant. Initial 
droplet drag (liquid acceleration) in supersonic flow tends 
to increase the static pressure (Table III). This regime 
is illustrated between taps No. 4 and #5 in run #32 and 
between taps #3 and #4 in run #44. 

In regime II evaporation is predominant. Evapora- 
tion in supersonic flow tends to decrease the static pres— 
sure (Table III). This regime is illustrated between tans 
#5 and #10 in run #32 and between taps #4 and #9 in run #44. 

In regime III the difference between water droplet 
temperature and gas stream temperature has decreased to the 
point where evaporation is no longer controlling, and wall 
friction becomes predominant, causing (Table III) an in- 
crease in static pressure. This is illustrated between 
taps #11 and #12 in run #32 and between taps #9 and #12 in 
run #44. 

The remaining portion of the curves in these two 
runs represent shock to subsonic flow with its accompanying 
pressure rise, followed by pressure increase induced by a- 


req increase in subsonic flow. 





The other runs in Figs. 7 and 8 represent operation 
with the shock in different locations, the shock patterns 
being clearly visible on the graphs. In runs #33 and #45 
the shock is in the diffuser, with the drop in static pres- 
sure between taps #12 and #13 being induced by the area in- 
crease in the diffuser. 

In run #26 (Fig. 7) the shock lies in the nozzle sec- 
tion. In this run regimes II and III take place in subgson- 
ic flow: regime I takes place in or near the shock gone. The 
static pressure rise between taps #4 and #10 is due to evan- 
oration: the drop between taps #10 and #12 is due to wall 


friction. 
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IX. Performance Comparisons 
A. Coefficient of Over-all Performance 
It will be found useful to write the fundamental e- 
quation of the A.T.P. (Section III) in "normalized" form 
by dividing through by the product of initial stagnation 


pressure and the square of the inlet Mach number: 


Po Pe h—) ($F dz 
Po. Mintet ' - 2m) fa To all D +7) 


In the search for the optimum inlet Mach number for 
the A.T.”. the basic consideration is the nroner balance 
between the contradictory requirements of high initial re- 
lative velocity and high initial temnerature difference 
between the air and water at the injection voint. Dowr- 
stream of the injection vnoint the local Mach number and 
therefore the gas stream temnerature may be controlled 
by area variation, as pointed out by Gavril (Ref. 7). 

Area variation used with moderation so as not to induce 
losses due to boundary layer phenomena has no effect upon 
stagnation pressure (Table TII). Controlled area variation 
in the design of a full-scale Aerothermo-ressor involves 
two contradictory requirements: high Mach number level 
through the duct (the factor M* in the fundamental 
equation), and high temperature difference between gas and 
water throughout the evaporative section. 


The answer to the question of balance between the 





requirements of high initial relative velocity and high 
initial temperature difference on the basis of small scale 
conatant-area experimental data lies in the magnitude of 
the last factor on the right-hand side of the funda- 
mental equation, rather than in the net over-all rise 
(or loss) in stagnation pressure, since the net over-all 
change in stagnation pressure can be seen to be strongly 
influenced by Mach number level throughout the length of 
the duct. 

By visualizing an integrated form of the normalized 
fundamental equation, it may be seen that for runs made 
in the same apparatus and having essentially the same ba 
profiles and M profiles (i.e., runs for which plots of 
Po/Po, vs. axial distance z along the duct would lie es- 
sentially along one and the same curve, and similarly for 


plots of M/Miniet vsS. Zz), the relative magnitudes of the 


f 


2 
quantity (Pop = Po, )/Po,Miniet or its equivalent (Pop/Po ma 


should yield the relative magnitudes inlet 
of the net effect of evaporation vs. friction over the duct 
length, as influenced by initial conditions. This quantity 
has in fact been commonly used on the A.T.P. Project as a 


coefficient of over-all performance. 


B. Supersonic Inlet vs. Subsonic Inlet 
Due to the inherent differences in the general flow 
character of the supersonic-inlet and subsonic-inlet cons- 


tant-area Aerothermopressors (opposite effects of evapo- 





ration and friction on Mach number, plus the presence of 
shock in one but not the other), comparison between a 
subsonic~inlet run and a supersonic~-inlet run by means 
of the over-all coefficient of performance defined above 
must be carried out with caution. 

In order to eliminate as many variables as vossible 
from the problem, comnarison of the sunersonic data was 
made with some vreviously unrevorted subsonic data re- 
corded by Mr. A.J. Erickson of the faculty at the Massa- 
chusetts Institute of Technology using the same test 
section and diffuser, rather than with the pvnublished data 
of Wadleigh (Ref. 4) which was taken on a 72" test section 
and 6° diffuser. 

Among Erickson's data are runs at 1500°R for in- 
let Mach numbers of 0.5 and 0.65 with different water-air 
ratios, utilizing Wadleigh's water injection nozzle No.3 
(Ref. 4) which injected water at the plane of the bell- 
mouth exit, His 1500°R run with best coefficient of over- 
all performance (run B-1-b-4 by his designation system) 
ie plotted in Fig. 10. 

As vreviously vointed out, run #26 (Fig. 7) re- 
presents suversonic-inlet with shock in the nozzle region 
followed by subsonic flow throughout the test section and 


diffuser, with consequent vn. and M profiles essentially 


O 
similar to Erickson's run in Fig. 10. For this reason 
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run #26 was chosen for the comparison. Its performance 
coefficient was comnuted utilizing a value of 1.47 for 
inlet Mach number, based on area ratio and measured 
pressure (Section VI). Performance coefficients for these 


two runs are tabulated below: 





Run Ms niet Por Po, Mi 

inlet 
MacKay #26 ey, , 868 —.061 
Erickson 
#B—1-b—4 0.65 .965 —- .083 


Sinee supersonic -inlet shows a substantially better 
coefficient of performance even while reflecting a shock 
loss, it appears that for an initial stagnation temmera- 
ture of 1500°, better evaporation rates are definitely 
obtainable at supersonic inlet Mach numbers. Certainly 
serious consideration of suversonic-inlet in the design 
of a large scale A.T. ?. seems warranted, especially with 
the possibility of using area variation to diffuse the 
supersonic flow to a lower Mach number before shock (hav-~ 
ing at the same time the effect of increasing the tempnera- 


ture differential between gas and water). 


C. Mach 1.5 vs. Mach 1.35 
In utilizing the performance coefficient of IX-A to 
compare Mach 1.5 with Mach 1.35, runs #27 and #39 were cho- 


sen, since these have shock regimes commencing at anvnroxi- 
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mately the same point in the cuct, and hence have somewhat 
Similar p, and M profiles. 

In the absence of dry nozzle pressure data for coni- 
cal tip No. 3, both performance coefficients for this compa- 
rison were computed using inlet Mach numbers corresponding 
to area ratio alone, as tabulated in Section VI. On this 


basis the performance data was as follows: 


(Pop/o,) — 1 





Run Minlet Dop/Po. 
Z M¢. 
inlet 
HOT 1.53 850 ~. 064 
739 1.34 .881 —. 066 


Since these two Mach numbers yield essentially equal 
performance coefficients while reflecting unequal shock losses, 
it seems reasonable to assume that the optimum inlet Mach 
number for an A.T.P. with controlled area variation may con- 
ceivably lie hizher than 1.47 (best estimate of flow Mach 
with tip No. 2) for inlet stagnation temperature of the 


order of 1500°R. 
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X. Suggestion for Further Study 

Because of the rather crude estimate of Mach numbers 
used in comparing Mach 1.5 with Mach 1.35 (nominal), no de- 
cision made on the basis of the computations in Section 
IX-C can be considered decisive. For a more accurate de- 
termination of the optimum supersonic inlet Mach number 
it would be advisable in any further small-scale testing 
to shorten the test section in order to obtain dry nozzle 
pressure data for inlet Mach numbers below 1.47. From runs 
#26. #33, and #45 (Figs. 7 and 8) it is evident that lit- 
tle would be lost by chopping off the present test section 
at tan No. 10, since at this point the undesirable influence 
of wall friction "takes over" from evaporation and becomes 
cont rolidine. 

Shortening the test section by this amount would. 
with the present water injectior nozgzgle_length, leave an 
effective constant area duct length of 24" extending from 
tan #4 (sunersonic nozzle exit) to tan #10. With an 
assumed f - .005 this gives 4fL/D - .226 and a value of 


M in the wéichorhood of 1.7. With this 


supersonic 
shorter test-section length and its corresvnonding reduction 
of minimum supversonic inlet Mach number for no friction 
choke from about 1.98 (Section IV) to about 1.7, it is 
possible (1) that dry supersonic nozzle data might be 


obtainable for conical tip No. 3 by auxiliary radial water 





mt 


injection downstream, and (2) a run for the Mach 1.35 
nozzle might be obtained having the same character as run 
#26 for the Mach 1.5 nogzle. With the annaratus used in 
this exneriment it was impossible with tip No. 3 to 

drive the shock any further upstream than run #38 (Fig. 8). 
since an increase in back pressure from this level re«# 


sulted in subsonic nozzle flow. 
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TABLE I 





LOCATION OF PERTINENT STATIONS ALONG THE DUCT 


Distance Downstream Distance From 


Station From Previous Bellmouth En 
Station (inches) trance (In. ) 

Stag. Pres. Probe #1 0 25% 

Bellmouth Entrance al O 

Nozzle Throat (Tap #2) 2 Ws 2148 

Tap #3 2 41/8 

A.T.P. Water Injection 1 578 5 3/74 

Nozzle Exit (tap #4) 3/8 6178 

Tap #5 4 10 178 

Auxiliary Water Injection 2 12 14% 

Tap #6 2 1418 

Tap #7 4 13° laze 

Tap #8 4 22 1/8 

Tap #9 4 26 IG 

Tap #10 \ 30 1/8 

Tap #11 4 34 1/78 

Tap #12 2 36 14 

Diffuser entrance 2 38.14 

Tap #13 27/8 M1 

Tap #14 8 4g 

Tap #15 8 Bui 

Tap #16 8 65 

Tao #17 8 73 

Tay #18 8 81 

Tap #19 8 89 

Tap #20 8 oF 

Tap #21 8 105 

Tap #22 8 ES 

Diffuser Exit 21/2 115 17 

Stag. Pres. Probe #23 9 Ley 172 

Quench tank entrance 19 143 1% 


* Indicates upstream 
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EVAPORATION, WALL FRIGTION, AND DROPLET DRAG 


Mach Number subsonic 
M 
supersonic 
Gas velocity 
V 
Pressure 
p ry 
Vall 
T 




























Area Evapora- Wall fric-+ Liquid accel- 






incresse |] tion pro- tion pro- | eration pro- 
produces duces duces duces 
(a) (b) (c) (a) 
















decrease | decrease(h)| increase | increase (e) 













increase | increase(h)| decrease | decrease (e) 





























subsonic decrease | decrease(h)} increase | increase(e) 








supereonic | increase | increase(h) | decrease | decrease(e) 
















subsonic increase | increase(h)| decrease | decrease(e) 










supersonic | decrease | decrease(h)| increase | increase(o) 





decrease(e) 





subsonic increase | decrease(h) | decrease 






supersonic | decrease | increase(h) |} increase | increase(e) 




















Gas Stagnation subsonic nil decrease nil decrease(f) 
Temperature 
supersonic decreasa nil decrease(f) 
Mixture Stag- subsonic decrease nil 
nation Temp. 
supersonic decrease nil 





Gas stagnation subsonic 


Pressure 





increase(h) | decrease | decrease(e) 





increase(h) decrease 


supersonic decrease(e) 






Mixture Stagne~ 
tion Pressure 


subsonic | Siersaeek A) decrease | decrease(z) 






incnseeSh) decrease 





decrease(g) 


gupersonic 


) Opposite effects for area decrease 

) Opposite effects for condensation 

) Opposite effects are impossible 

) When y<1, d¥g>0O; when y>1l, dVg<0O (See Ref. 2) 

) Dependent upon magnitude of y for liquid decelera- 
tion (See Ref. 2) 

) Opposite effect for liquid deceleration 

) Same effect for liquid deceleration 

) Based on # only, and generally correct for A 
in excess Of two; otherwise effects are indeter= 
minate (See Ref. 2) 
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